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We study properties of the Abelian monopoles in the Maximal Abelian projection of the three 
dimensional pure SU(2) gauge model. We match the lattice monopole dynamics with the continuum 
Coulomb gas model using a method of blocking from continuum. We obtain the Debye screening 
length and the monopole density in continuum using numerical results for the lattice density of 
the (squared) monopole charges and for the monopole action. The monopoles treated within our 
blocking method provide about 75% contribution to the non-Abelian Debye screening length. We 
also find that monopoles form a Coulomb plasma which is not dilute. 

PACS numbers: 11.15.Ha,14.80.Hv,11.10.Wx 



I. INTRODUCTION 



According to the dual superconductor mechanism the confinement of quarks in non-Abelian gauge theories is 
caused by the monopole-like configurations of the gluonic fields. These configurations can be identified with the help 
of the Abelian projection method 0. The basic idea behind this method is to fix partially the non-Abelian gauge 
symmetry up to an Abelian subgroup. If the original non-Abelian gauge group is compact (like the SU(N) group) 
then the residual Abelian symmetry group is compact as well. The compactness of the Abelian group guarantees the 
existence of the monopoles. 

In the low temperature phase of the four dimensional SU(N) gauge model the monopoles are condensed in 
agreement with the dual superconductor mechanism The condensation leads to the appearance of the dual 
Meissner effect and, as a result, to the formation of the chromoelectric string between fundamental sources of the 
chromoelectric field (quarks). Consequently, the quarks get confined. Moreover, the Abelian monopoles are in the so- 
called Maximal Abelianprojection Q make to make a dominant contribution to the zero temperature string tension^ 
(for a review, see Ref. ||). 

At the critical temperature, T — T c , the monopole condensate disappears and at higher temperatures the quarks 
are no more confined. The vacuum in the deconfincmcnt phase is filled by the static monopoles which can not lead to 
Oh' the confinement of the static quarks. However, the absence of the color confinement at high temperatures does not 
mean that the high temperature physics is perturbative. Indeed, even at T > T c the quarks running in the spatial 
directions are still confined due to existence of the " spatial string tension" (this is a coefficient in front of the area 
term of large spatial Wilson loops) . The spatial string tension is a non-perturbative quantity which is dominated by 
J^j ' contributions of the static monopoles according 1 to Ref. . 

The physics of the static monopoles in the high temperature SU(2) gauge model was investigated in Ref. Q using the 
method called blocking from continuum (BFC). This method resembles the idea of the blocking of the continuum fields 
to the lattice ^(j. I n general the BFC procedure allows us to match lattice results to the continuum model without 
going to a deep continuum limit. For example, the BFC method clearly shows 9] that in the continuum limit the static 
monopoles in the high temperature SU(2) gluodynamics are described by the three dimensional Coulomb gas model. 
Another example is the zero-temperature 4D SU(2) gluodynamics in which the value of the monopole condensate can 
be obtained with the help of the BFC method [ 1 1| . Being encouraged by these results, we apply the BFC method to 
the monopoles in the three-dimensional SU(2) gluodynamics. The confinement mechanism in this model based on the 
Abelian monopole dynamics was previously discussed both within analytical [T^ and numerical |!3lll4j frameworks. 

The plan of the paper is as follows. In Section [H] we briefly recall the results of the Ref. [jj for the lattice density 
of the (squared) monopole charge and for the monopole action expressed via parameters of the continuum Coulomb 
gas model. In Section ITTT1 we study the lattice action and the lattice density of the monopoles in the pure 3D SU(2) 
gauge model. Using the BFC method we get density of the monopoles and the monopole contribution to the magnetic 
screening length in the continuum limit. Our conclusion is presented in the last Section. 



1 In a different scenario the "spatial confinement" problem is suggested to be caused by magnetic thermal quasi-particles. It seems 
plausible that these magnetic excitations are correlated with the Abelian monopoles. 
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II. LATTICE MONOPOLES FROM CONTINUUM MONOPOLES 



The standard way to identify the lattice monopoles in Monte Carlo simulations is to use the DeGrand-Toussaint 
construction which calculates the magnetic flux coming out of lattice 3D cells (cubes). The magnetic charges 
obtained in this way are conserved and quantized. The properties of the lattice monopoles should obviously depend 
on the physical size, b, of the lattice cell (below we call these lattice objects as "lattice monopoles of the size b"). To 
get the properties of the monopoles in continuum one should send the size of the lattice cells to zero, b — ► 0, what is 
usually a difficult numerical problem. The BFC method allows to get the properties of the monopoles in continuum 
(" continuum monopoles" ) using the results obtained on the lattice with a finite lattice spacing. 

The idea behind the BFC method is to treat each lattice 3D cell as a "detector" of the magnetic charges of the 
continuum monopoles. If the continuum monopole is located inside a lattice 3D cell then the DeGrand-Toussaint 
method detects an existence of a lattice monopole inside this cell. If the size, b, of the lattice cell is finite then two 
or more continuum monopoles may be located inside the cell. The fluctuations of the monopole charges of the lattice 
cells must depend on the properties of the continuum monopoles. As a result, the lattice observables - such as the 
vacuum expectation value of the lattice monopole density (or, the lattice action) - must carry information about 
dynamics of the continuum monopoles. The observables should depend not only on the size of the lattice cell, b, but 
also on features of the continuum model which describes the monopole dynamics. 

To avoid misunderstanding we would like to stress from the very beginning the difference between various lattice 
monopole sizes b. As we mentioned above, we call the size of the lattice 3D cell - used for detection of the monopoles 
- as "the size of the lattice monopole". This should be distinguished from the physical radius, ro, of the monopole 
core which is obviously independent on the size of the lattice " detector" . In this paper we disregard the existence 
of the monopole core and consider the continuum monopoles as point-like objects. Thus, we get the finite-sized 
lattice monopoles by blocking the point -like continuum monopoles to the lattice. 

Below we describe briefly the techniques which appear in the BFC method. Let us consider a three-dimensional 
lattice with a finite lattice spacing b which is embedded in the continuum space-time. The cells of the lattice are 
defined as follows: 

C, = \b(a i -l)<x i <b(a i +l), i= 1,2,3}, (1) 



2/ ~ V 2 

where Sj is the lattice dimensionless coordinate and corresponds to the continuum coordinate. 
The magnetic charge, k a , inside the lattice cell C s is 

k s = f d 3 x p{x) , p(x) = ]>> <5 (a) [x - x (a) ) , (2) 

where p{x) is the density of the continuum monopoles, x a and q a is the position and the charge (in units of a 
fundamental magnetic charge, gu) of a th continuum monopole. In three dimensions the monopoles are instanton-like 
objects and the monopole trajectories have zero dimensionality (points). Basic properties of the lattice monopoles 
are the same as the ones of the continuum monopoles: the charge k s is quantized, k s € "Z, ' , and conserved in the 



three-dimensional sense: 

^2k s = d 3 x p(x) = . 
seA I 

Here A and V denote the lattice and the continuum volume occupied by the lattice, respectively. 

In the BFC method various properties of the lattice monopoles such as the lattice monopole density, correlators of 
the lattice monopole charges, the effective lattice monopole action etc. can be calculated in a continuum monopole 
model using Eq. @ as a definition of the lattice monopole charge. In this paper we suppose that the dynamics of the 
continuum monopoles is governed by the 3D Coulomb gas model: 



N=0 



N . "1 2 N 

II / d 8 *<-> £ exp{-^ £ q a a b D{x^ - *«)} . (3) 



9a=±l. 



.,6—1 

a=fb 



The Coulomb interaction in Eq.@ is represented by the inverse Laplacian D, —dfD(x) = S^ 3 '(x), and £ is the fugacity 
parameter. 
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The model @ does not exist without properly defined ultraviolet cut-off. Indeed, the self-energy of the point-like 
monopoles is a linearly divergent function. As a result, the fugacity must be renormalized, £ rcn = C ■ ex p{ffm/ (87r r )}, 
where ro is the ultraviolet cut-off. In our case this cut-off is given by the size of the monopole core which is of the 
order of 0.05 fm at zero temperature [16j. For simplicity we omit the subscript "ren" in the renormalized fugacity 
below. 

The magnetic charges in the Coulomb gas (01 are screened: at large distances the two-point charge correlation 
function is exponentially suppressed, (p(x)p(y)) ~ exp{ — \x — y\/Xo}- Here Ad is the Debye screening length |17| . 



(4) 



which is inversely proportional to the Deb ye s creening mass, Mo = \ D . The density of the continuum monopoles in 
the leading order is related to fugacity as 17] p = 2£. 

In Ref. [9j the lattice monopole action, S mon (k), and the v.e.v. of the squared magnetic charge, (fc 2 ), were calculated 
starting from the Coulomb gas model ©. In the low-density approximation the leading order contribution to the 
(squared) density of the monopole charges is Q : 



( fc2 ( fe )) = T3 (E fc2 ( fi ))^ / d " x [ d 3 y(p(x)p(y))=pb 3 P(M D b), 



where in the thermodynamic limit (an infinite- volume lattice) the function P is 



P(p) = l-p 2 



1 



(2tt) 3 q 2 + M 1 



n 



2sin( ft /2) 



-i 2 



<1< 



The finite-volume analogue of Eq. @ can be obtained by the standard substitution: 

oo 

271-fci f da; 1 \ - 

77 



q, 



2tt 



(5) 



(6) 



(7) 



ki£2 



where L; is the lattice size (in units of the lattice spacing) in i th direction. 

The reason why the density of the squared magnetic charges, (fc 2 ), is used instead of the standard definition of the 
density, (|fc s |), is simple: an analytical treatment of (\k s \) is obviously much more difficult compared to that of the 
quantity (fc 2 ). 

The explicit behavior of the lattice monopole density (0 as the function of b can be obtained in the limit of large 
lattice monopoles 0, 



6> Ad, 



(8) 



6< A D , 



(9) 



(k 2 ) = C lP \ D b 2 - \ l + 0((\ D /b) 
as well as in the case of small monopoles, 

(k 2 ) =pb 3 - [l + C 2 p(b/X D ) 2 +o((b/\ D )' 

where C x « 2.94 and C* 2 « 0.148. 

The proportionality of the density (fc 2 ) to b 2 in large-6 region has a simple explanation 9]. In a random gas of 
continuum monopoles we would obviously get (fc 2 ) ~ pb 3 . Due to the Debye screening in the Coulomb monopole gas 
the monopoles separated from the boundary of the cell by the distance larger than Ad, do not contribute to (fc 2 ). 
Consequently, the b 3 proportionality for the random gas turns into Xd^ 2 in the Coulomb gas and we get (fc 2 ) cx p Xob 2 . 

In the small b region the density of the squared lattice monopole charges is equal to the density of the continuum 
monopoles times the volume of the cell. This is natural, since the smaller volume of the lattice cell, b 3 , the smaller 
chance for two continuum monopoles to be located at the same cell. Therefore each cell predominantly contains not 
more that one continuum monopole, which leads to the relation fc 2 = \k s \ = 0, 1. As a result we get (fc 2 ) — > pi a tt(b) — > 
pb 3 in the limit b — > 0. 

Analogously one can get that in the large-6 region the action of the monopole is given by the lattice Coulomb 
action 



*5 'mon (fc) 



pXi 



1 

ft 2 



s,s' 



k s D s 8 > k s i + 



b > Ad 



(10) 



where the coefficient in front of the lattice monopole action is inversely proportional to b 2 . 
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III. NUMERICAL RESULTS 

We have simulated the pure SU(2) gauge model with the standard Wilson action S = — l/2^ p Trf/p, where 
Up is the plaquette matrix constructed from the gauge link fields, U\. We have generated 200 configurations of the 
gauge fields for each chosen value of the coupling constant, f3 = 2.083, 2.5, 3, 3.47, 3.75, 4.5, 5, 6, 6.56, 9, 12, 14.5, on the 
lattice 48 3 . To study the Abelian monopole dynamics we perform Abelian projection in the Maximally Abelian (MA) 
gauge -4] for each SU(2) configuration. The MA gauge fixing condition is the maximization of the quantity R, 

wgKRpM], R[U} = TvJ2[U^s)a 3 Ul(s + fi)a 3 }. (11) 

under the SU{2) gauge transformations, U — > = &UQ. The gauge fixing condition 1|11|) is invariant under 

an Abelian subgroup of the group of the SU(2) gauge transformations. Thus the condition (|llfl corresponds to the 
partial gauge fixing, SU{2) — » U(l). 

After the MA gauge fixing, the Abelian, {m p (s)}, and non-Abelian {U^s)} link fields are separated: 

%{s) = C,( sMs) , CM = ( VT ^f W ^T^W ) • - ( e-i« ) • < 12 > 

The vector fields C M (s) and u M (s) transform like a charged matter and, respectively, a gauge field under the residual 
U(l) symmetry. Next we define a lattice monopole current (DeGrand-Toussaint monopole) Abelian plaquette 

variables 9 )JlV {s) are written as 

M«) = 0m(*) + 0*(* + + *)-*«'(*)> (-4tt < 6^{s) < 4tt). (13) 

It is decomposed into two terms using integer variables ro^^s): 

9 flv (s) = e f , v (s) + 2Trn flv (s), (-tt < 6^(s) < tt). (14) 

Here 6^ u {s) is interpreted as an electromagnetic flux through the plaquette and n M „(s) corresponds to the number of 
Dirac string piercing the plaquette. The lattice monopole current is defined as 

k ( s ) = \^v P ad v n prj {s + /}) . (15) 

In order to get the lattice density for the monopoles of various sizes, b, we perform numerically the blockspin 
transformations for the lattice monopole charges. The original model is defined on the fine lattice with the lattice 
spacing a and after the blockspin transformation, the renormalized lattice spacing becomes b = na, where n is the 
number of steps of the blockspin transformations. The continuum limit is taken as the limit a — > and n — > oo for a 
fixed physical scale b. 

The monopoles on the renormalized lattices ("extended monopoles", Ref. [l^) have the physical size 6 3 . The charge 
of the n-blocked monopole is equal to the sum of the charges of the elementary lattice monopoles inside the n 3 lattice 
cell: 

k ( - n) (s)= ^2 k(ns + ifi + jp + lp). 

i,j,t=0 

For the sake of simplicity we omit below the superscript (n) while referring to the blocked currents. We perform the 
lattice blocking with the factors n = 1 ... 12. All dimensional quantities below are measured in units of the string 
tension, a, the values of which are taken from Ref. [I^.l20| . 

We show the density of the squared monopole charges (normalized by the factor b~) in Figure ^ as a function of the 
scale b for various blocking factors, n. One can easily see from this Figure that the density shows only an approximate 
scaling: the data depends not only on b = an, but also on the fine lattice spacing, a, and the blocking factor, n, 
separately. The violation of the 6-scaling - being much bigger then the error bars of the data - strongly affects the 
physical results. 

The origin of the 6-scaling violation is in the presence of the ultraviolet lattice artifacts which express themselves 
in the form of the tightly bound monopole-anti- monopole pairs (magnetic dipoles). These dipoles are living on the 
fine lattice and typical distance between the constituents of the these dipoles is of the order of the fine lattice spacing, 
a. Thus, in the continuum limit these dipoles must disappear. However, on the finite lattices they affect the results 
drastically. 
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FIG. 1: The density of the squared lattice monopole charges, (kg), divided by b 2 as a function of 6 for various blocking steps n 
(b is given in units of the string tension, a). 

In order to get rid of the ultraviolet artifacts we have removed the tightly-bound dipolc pairs from all configurations 
using a simple numerical algorithm. Namely, we remove a magnetic dipole if it is made of a monopole and an anti- 
monopole which are touching each other (i.e., this means that the centers of the corresponding cubes are located at 
the distance smaller or equal than \/3a). Note, that we first apply this procedure to the elementary a 3 -monopoles, 
and only then we perform the blockspin transformations. Below we discuss the results obtained for the monopole 
ensembles with the artificial UV-dipoles removed. 
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FIG. 2: The density of the squared monopole charges, (fef), with the UV dipoles removed. The density is normalized (a) by b 2 
and (b) by b 3 . The fits by the function JSJ are shown by dashed lines for each value of the blocking step, n. 

In Figure Ufa) we show the density of the squared monopole charges as a function of the scale b for various values 
of the blocking factor, n. The density is normalized by the factor b 2 . One can see that the 6-scaling violations are 
very small. As the blocking size b increases the slope of the ratio (k 2 )/b 2 becomes less steep. This behavior is in a 
qualitative agreement with the prediction from the Coulomb gas model which says that in the high-6 limit the 
ratio (k 2 )/b 2 should converge to a constant. 

On the other hand the quantity {kg) should be proportional to b 3 in the small-6 region JjJ}. One can indeed observe 
from Figure|3b) that the ratio (k 2 )/b 3 does tend to a constant at small b. Note that there is a small scaling violation 
in this region which is due to the presence of the lattice artifacts at the scale b ~ a. In order to get artifact-free 
results we will use below large-6 monopoles. 
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The values of the parameters of the Coulomb gas model in the continuum limit, Eq. can be obtained by fitting 
the numerical results for (k 2 ) by the theoretical prediction JSJ),J7|. Technically, for each value of the blocking step, n, 
we have a set of the data corresponding to different values of the lattice coupling /3, and, consequently, to different 
values of b = n ■ a{f3). Note that by fixing n we simultaneously fix the extension of the coarse lattice, L/n, in units of 
b. The size of the coarse lattice enters in Eq. (JJJ. We fit the set of the data for the fixed blocking step n. The best 
fit curves are shown in FiguresE^a) and (b) by dashed lines. The quality of the fit is very good, \ 2 /d.o.f. ~ 1 . 

The density fits give the values of the continuum monopole density, p^ , and the Debye mass, Mjy , which are 
shown in Figures E3a,b) as a function of blocking step n. All these results are given in units of the string tension. 
The influence of the finite lattice spacing should vanish in the limit of large 6, or, in our case, in the limit of large n: 




12 n 




12 n 



FIG. 3: (a) The density of the continuum monopoles, p, and (b) the Debye screening mass, Md, obtained with the help of the 
fits of the n-blocked squared monopole density by function JSj . The large-n extrapolation 1161 1 is shown by solid lines. 

lmin^oo £K n ) = O ph , where O stands for either p or M D and the superscript "ph" indicate the artifact-free physical 
value. We found that for n > 2 the dependence of both p and Mp on the blocking size n can be approximately 
described as 



(n) = ph + const . n -a 



(16) 



The extrapolation (|16[) is shown in Figure [21 We get the physical values for the monopole density p and the Debye 
screening mass Mn coming from the Coulomb gas model (here and below we omit the superscript "ph" for the 
extrapolated values): 



p / a i/z = 0.174(2) , M D /a i/2 = 1.77(4) . 

The value of Md may be treated as the "monopole contribution to the Debye screening mass" 
The self-consistency check of our approach can be done with the help of the quantity 



C* = 



M D a 



(17) 



(18) 



which is known to be equal to eight (C CG = 8) in the low density limit of the Coulomb gas model 01- I n Figure 0] 
we plot our numerical result for C as a function of n. 
The large-n extrapolation (|16(l gives 



C= 10.1(1). 



c 

CCG 



1.26(3). 



(19) 



The quantity C is about 25% larger than the one predicted by the Coulomb gas model in the low monopole density 
approximation, C^, G = 8. The discrepancy is most likely explained by the invalidity of the assumption that the 
monopole density is low. Indeed, the low-density approach requires for the monopole density to be much lower than 
a natural scale for the density, g 6 (remember that the coupling g has the dimensionality mass 1 / 2 ). The requirement 
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FIG. 4: The same as in Figure|!|]but for the ratio 1181 . The dashed line corresponds to the low density limit of Coulomb Gas 
model E3, C CG = 8. 

p <C g 6 can equivalently be reformulated as p/Mp ^> 1, which means that the number of the monopoles in a 
unit Debye volume, A|> = M£, 3 , must be high. Taking the numerical values for p and Mp from Eq. (|17|l we get: 
p/Mjj w 0.03 <C 1. Thus, the low-density assumption is not valid in the 3D SU(2) gluodynamics. However, the 
discrepancy of 25% observed in the quantity C, Eq. H19[). is a good signal that the Coulomb gas model may still 
provide us with the predictions valid up to the specified accuracy. 

One can compare our result for the monopole density, Eq. (|17H . with the result obtained by Bornyakov and Grigorev 
in Ref. 0, p BG = 2~ 7 (1 ± 0.02) g 6 . Using the result of Ref. Jo = 0.3353(18) g 2 , we get the value p BG /cr 3 / 2 = 
0.207(5), which is close to our independent estimation in the continuum limit l|17fl : p/p BG = 0.83(4). The result of 
Ref. [l3| is about 20% higher than our estimation for the monopole density. Thus, although the condition of the low 
monopole density approximation is strongly violated, the BFC method (based on the dilute gas approximation) gives 
the value of the monopole density which is consistent with other measurements. 

It is interesting to compare the result for the screening mass (|17f) with the lightest glueball mass measured in 
Refs. 0,|2(j, Mq++ = 4.72(4) Jo. In the Abelian picture, the mass of the ground state glueball obtained with the 
help of the correlator, 

KM F ap( R )) = COnSt - e ~ M ° ++ R + ■ ■ ■ > 

must be twice bigger than the Debye screening mass, 2Md/Mq++ = 1, where the Debye mass is given by the following 
correlator 

(i^(Q) F^(R)) = const. e' M ° R + . . . . 

The comparison of our result with the result of Refs. [H |2(J gives 2Md/M ++ = 0.75(4). The deviation is of 
the order of 25% similarly to case of the quantity C. 

Let us also compare our result for the monopole contribution to the Debye screening mass, Eq. Ijl7|l . with the direct 
measurement of the Debye mass in 3D SU(2) gauge model made in Ref. [22, ^d^ 2 /V& = 1-39(9). The values agree 
with each other within the 25 per cent: mo/m^ 1 ^ — 1.27(11). Approximately the same accuracy is observed in the 
four-dimensional SU(2) gauge theory for the monopole contribution to the fundamental string tension [22|. 

Finally, we have made the cross-check of our result by a numerical calculation of the effective monopole action. We 
used an inverse Monte-Carlo algorithm described in Ref. [23] and we used the following form of the trial action: 

S m on{k) =fl-) J kl + fcoul • k s D s , s , k s > . (20) 
s s,s' 

The choice of this type of the trial action is motivated by the following reasons. From the point of view of reliability 
of the results the data is most valuable in the large-6 region since the subtraction of the ultraviolet monopoles can 
not change the infrared physics. In this region the action is expected to be given by the Coulomb term l|10|) . On the 
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other hand the subtraction of the ultraviolet monopoles must affect the ultraviolet (or, local) terms in the monopole 
action. Thus we added to the trial action H20|) the most local term with the coupling fx. The role of this term is to 
take into account the effects which are caused by the monopole subtraction. 

We depict the coupling /i and the product b 2 fcoul as a function of b in Figures |SJa,b), respectively. The strong 
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FIG. 5: The coupling /i and the product b 2 fcoui of the action 121 H vs. b for each value of the blocking step, n. 

coupling region with (3 $S 3.0 was not included into these figures due to large finite size effects. Indeed, in this case 
a^fa > 0.6 which implies that aMp « 1 (here we used the value of Md = A^ 1 obtained from the fits of the monopole 
density l(T7Jll. 

Note that small-n results are more sensitive to the subtraction of the ultraviolet magnetic dipoles. This is the 
reason for noticeable ^-scaling violations around b ~ 0.5a- 1 / 2 for n = 1 in both Figures EJa,b). 

According to the analytical prediction l|10|l the product of the coupling fcoui and b 2 should tend to a constant in 
the limit b^$> Xu while the /i-term should disappear from the action. Numerically, the coupling fi is diminishing as 
b increases, as seen in Figure EJa). One can also check numerically that this coupling is vanishing quicker than 1/b 2 
at b > 1.5er~ 1//2 . So in the large-6 region we are indeed left only with the lattice Coulomb term in the action. If the 
scale b is expressed in units of the string tension then the product b 2 fcoui{b) tend to the quantity fl^jl . Figure OJb) 
shows that large-& plateau in distribution of b 2 fcoui(b) starts from b > 2er~ 1 / 2 . Averaging the available data with 
b ^ 2a- 1 / 2 we get C = 10.2(3) which is in an excellent agreement with the result l|19|) obtained from the monopole 
action. 



IV. CONCLUSIONS 



We have shown that the dynamics of the Abelian monopoles in the three-dimensional SU(2) gauge model can be 
described by the Coulomb gas model. Using a novel method, called the blocking of the monopoles from continuum, we 
calculated the monopole density and the Debye screening mass in continuum, Eq. I|17|) . using the numerical results for 
the (squared) monopole charge density. The self-consistency of the results was checked by the independent analysis 
of the lattice monopole action. We conclude that the Abelian monopole gas in the 3D SU(2) gluodynamics is not 
dilute. Nevertheless, the continuum values of the monopole density (p = 0.174(2) a 3 / 2 ) and the Debye screening mass 
{Mjj = 1.77(4) a 1 / 2 ) - obtained with the help of the dilute monopole gas model - are consistent within the accuracy 
of 25% with the known data obtained from independent measurements. 
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